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Abstract 

In this paper we describe in full details a new family of recently found exact solutions of rela- 
tivistic, perfect fluid dynamics. With an ansatz, which generalizes the well-known Hwa-Bjorken 
solution, we obtain a wide class of new exact, explicit and simple solutions, which have a remark- 
able advantage as compared to presently known exact and explicit solutions: they do not lack 
acceleration. They can be utilized for the description of the evolution of the matter created in high 
energy heavy ion collisions. Because these solutions are accelerating, they provide a more realistic 
picture than the well-known Hwa-Bjorken solution, and give more insight into the dynamics of 
the matter. We exploit this by giving an advanced simple estimation of the initial energy density 
of the produced matter in high energy collisions, which takes acceleration effects (i.e. the work 
done by the pressure and the modified change of the volume elements) into account. We also give 
an advanced estimation of the life-time of the reaction. Our new solutions can also be used to 
test numerical hydrodynamical codes reliably. In the end, we also give an exact, 1-1-1 dimensional, 
relativistic hydrodynamical solution, where the initial pressure and velocity profile is arbitrary, and 
we show that this general solution is stable for perturbations. 

PACS numbers: 24.10.Nz,47.15.Hg 
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INTRODUCTION 



Hydrodynamics is in a unique position among the different branches of theoretical physics: 
it rehes only on the assumption of local charge and energy-momentum conservation, and on 
the concept of local thermal equilibrium. The hydrodynamical equations do not have an 
internal scale, in contrast to other theories. Consequently, among the possible applications 
of hydrodynamics we find the largest known physical systems, systems of intermediate size 
and the smallest and shortest events of nature as well. For example, it is well known, that 
the evolution of the whole Universe (the Hubble flow) can be treated in a hydrodynamical 
manner. The applications in the ordinary sizes and temporal extents are known since the 
birth of hydrodynamics. Considering the smallest temporal and spatial scales: as a surprising 
new result, it has been understood recently, that the new form of nuclear matter created in 
high energy heavy ion collisions can be considered as an (almost) perfect fluid — because 
its properties (like the single-particle spectra, elliptic flow data and two-particle correlation 
functions) are best explained in terms of models based on perfect fluid hydrodynamics. 
From this point of view, we can say that hydrodynamics is simple and efficient, and it has 
applications from the very small to the very large scales. 

The hydrodynamical equations, however, are highly nonlinear even in the nonrelativistic 
case, this is the reason for the appearance of chaos, instability and other beautiful flow pat- 
terns in realistic flows. Recent interest in hydrodynamics emerged mostly from the challenge 
to explain the properties of the nuclear matter in high energy collisions, it is clear that in 
this context one has to deal with not only hydrodynamics, but relativistic hydrodynamics. 
Relativity makes the case even worse: relativistic terms (such as speed-addition and proper- 
time factors) make the equations much more complicated as compared to the nonrelativistic 
case. From this point of view, hydrodynamics is a horrendously complicated area of physics, 
which can be seen also from the intensity of recent efforts — both analytical and numerical 
— to find relevant solutions for the relativistic hydrodynamical equations. 

There are only a few exact solutions for these equations. One (and historically the first) 
is the famous Landau-Khalatnikov solution discovered more than 50 years ago [1-3]. This 
is a 1-1-1 dimensional solution, and has realistic properties: it describes a 1-1-1 dimensional 
expansion, and does not lack acceleration. Landau was able to calculate approximately an 
important observable, the rapidity distribution of the flowing particles from this solution. 
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and it was found to be an (approximately) Gaussian distribution. Recent measurements at 
the RHIC particle accelerator (BNL, USA) suggest that this prediction is not at big vari- 
ance with the observations made even at the highest presently attainable collision energies. 
However, considering its formulation, the Landau-Khalatnikov solution is a very unpleasant 
one: it is an implicit solution, the independent variables (the space and and the time co- 
ordinate) are given by extremely complicated integral formulas, involving the fluid rapidity 
and temperature. 

Another renowned relativistic hydrodynamical solution is the Hwa-Bjorken solution [4, 5], 
which is a simple, explicit and exact, but accelerationless solution. This solution (in its 
original form, where no conserved baryonic charges are present) is boost-invariant, and from 
this it follows that the rapidity distribution should be constant. This prediction would be 
valid for the limiting case of infinitely high center of mass collision energy. Although such 
a fiat ,, plateaux" (a narrow domain of constant rapidity distribution around mid-rapidity) 
is indeed seen in recent rapidity measurements at RHIC, the observed distribution is far 
from constant, so Bjorken's approximation fails to describe the data. However, the solution 
allowed Bjorken to obtain a simple estimate of the initial energy density reached in high 
energy reactions from final state hadronic observables. (This is why mostly only Bjorken is 
credited with this solution, though Hwa discovered it independently almost 10 years before.) 
The Bjorken-estimation is widely used, but it is well known, that the acceleration of the 
matter influences the energy density estimation. In order to obtain more realistic estimation, 
one should use explicit and exact, but accelerating solutions of relativistic fluid dynamics. 

In this paper we present a detailed description of such new simple, exact and explicit 
solutions of relativistic fluid dynamics. The essential features of these solutions were first 
presented in two short Letters in refs. [6, 7]. To our best knowledge, these solutions were 
the first explicit solutions with relativistic acceleration, following the implicit and accel- 
erating Landau-Khalatnikov solution. (There are other interesting solutions to relativistic 
hydrodynamics, see e.g. refs. [8-11], but these solutions are all accelerationless.) 

Apart from the purely theoretical point of view of finding new exact solutions for nonlinear 
differential equations, our new solutions are important due to other reasons, too. First, they 
provide an advanced estimate of the initial energy density and life-time of high energy 
reactions, we focus on these applications in this paper. Second, they can be applied to 
test numerical solutions of relativistic hydrodynamics: no finite, accelerating relativistic 
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solutions were available before for 1+3 dimensional tests. 

The organization of the paper is as follows. We discuss the equations of relativistic 
hydrodynamics, then present the newly found solutions. We then calculate the rapidity 
distribution, and use it to estimate the initial energy density and the life-time. In the four 
Appendices, details of calculations are summarized. In Appendix A we derive the solutions 
themselves, and in Appendix B we prove their uniqueness (in some sense discussed there). 
Appendix C presents the calculation of the rapidity distribution, and finally, in Appendix 
D we present a general solution and investigate the stability of the solutions for a specific 
choice of the equation of state. 

NOTATION AND BASIC EQUATIONS 

In this section we specify the notations used throughout the paper, and derive the equa- 
tions of relativistic hydrodynamics. 

Notation 

We use three-dimensional as well as four-dimensional notations depending on which one 
is more convenient. In four-dimensional notation the independent variables are the com- 
ponents of the x'^ coordinate four- vector: x'^ = {t,r), with t being the time and r the 
coordinate three-vector, r = {rxify^Tz). The metric tensor is denoted by g^^, we use the 
g^u = diag {1, —1, —1, —1) sign convention. The four-velocity field is denoted by m^, nor- 
malized to unity, m^-u^ = 1, that is, we treat the speed of light as 1. The three- velocity v is 
defined as = 7 (1, v), with 7 = (1 — v^)^^^'^. We denote the so-called comoving derivative 
bv -• = + vV 

The thermodynamical quantities are the following: e is the energy density (including 
rest energy contribution), p is the pressure, w = e + p is the enthalpy density, T is the 
temperature, a is the entropy density. When there are some (conserved or non-conserved) 
charges present, we denote them by rij, and the corresponding chemical potentials by /ij. 
All the above densities are defined in the local rest frame, so these densities are scalar 
quantities, while for example the entropy density in the laboratory frame is defined as the 
0-th component of the entropy current vector, au^. 
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The equations of relativistic hydrodynamics 



Assuming perfect fluid (tfiat is, neglecting heat conductivity, bulk and shear viscosity) the 
equations of relativistic hydrodynamics are obtained by Landau's heuristic argumentation: 
we guess the form of the energy-momentum tensor T^j, in the local rest frame: T^^ = 
diag {e,p,p,p), and from this we have in arbitrary frame 

T^^ = wu^u^ - pg^y. (1) 

The equations then follow from the energy-momentum conservation law d^T^^ = 0. Substi- 
tuting eq. (1) and projecting orthogonal and parallel to u'*, we obtain the relativistic Euler 
equation and the energy conservation equation as 

wu^d^u^" = {g^^ - u^'uP) dpp, (2) 
wd^u^" = -u^'dpE. (3) 

The general form of the charge conservation equations (for one type of charge) is as follows: 

J]/i,9^(niM'^) = 0, (4) 

i 

this is valid for the case of non-conserved charges, as their chemical potentials vanish, for 
one conserved charge n with non-vanishing chemical potential, as the conservation law is 

{nu^) = 0, and also for various mixtures. (For instance, in case of baryonic or electric 
charge, particles and antiparticles carry opposite charges, and they chemical potentials are 
the same but of opposite sign, so eq. (4) is valid also for this case.) 

For latter usage we write down these equations also in a three-dimensional notation. The 
relativistic Euler equation, the energy conservation equation and the continuity equation 
(for one conserved charge) are 

w dv dp\ 

'Vp + v-^ , (5) 



dt V 

The thermodynamical quantities obey general rules, which follow directly from thermody- 
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namics: 

de = Tdfj + /ijdrij, (8) 

i 

W = Tcr + ^ fXiUi, (9) 

i 

dp = adT + njd^j. (10) 

i 

The last equation is a consequence of the first two ones. From eqs. (8) and (10) and eq. (3) 
one can derive the entropy conservation equation: 

d, (an^) = 0. (11) 

This is a consequence of the special choice of the energy-momentum tensor. The perfectness 
of the fluid can be defined alternatively by requiring that eq. (11) holds. 

In order to get a closed set of equations, we have to specify the Equation of State (EoS) 
of the flowing matter. In the following we shall investigate a case, which resembles to that 
of ideal gas, that is 

e = up. (12) 

The K constant is 3 for a three-dimensional, ultra-relativistic ideal gas, but now we retain the 
more general possibility of arbitrary n. With the choice of eq. (12) the charge conservation 
equation decouples from the Euler and the energy conservation equations. That is, in order 
to find new solutions, one has to deal with v and p only, because eqs. (5) and (6) contains 
only these variables, and charge densities occur in the calculation only after one found proper 
solutions for p (t, r) and v (t, r). 

An important property of this specific EoS is, that if one finds a solution with this EoS, 
then these solutions can be generalized for a bag-model type of EoS in a straight-forward 
manner, as the following replacements leave the hydro dynamical equations invariant: 

e ^ e + B, (13) 
p p - B. (14) 

The value of the bag constant, B cancels from eqs. (5) and (6), that is, one can introduce 
an arbitrary value of the bag constant. This does not change the time evolution of the 
entropy density or the flow field as long as it is a constant, independent of position and 
time. Because this possibility is simple, we do not explicitly write down the bag constant in 
the following calculations, but emphasize, that this can be done without difficulty. 
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Rindler coordinates 



We discuss below exact solutions of relativistic perfect fluid hydrodynamics in 1+1 di- 
mensions and spherical solutions in 1+d dimensions as well (the 1+1 dimensional case is 
a special case of this latter one). The number of spatial dimensions is denoted by d. The 
notation r stands for the spatial coordinate in 1 + 1 dimensions, and for the radial coor- 
dinate in 1+d dimensions. All quantities depend only on t and r, and the velocity field is 
radial. We use the well-known Rindler coordinates r and rj as independent variables, which 
naturally fit to the Hwa-Bjorken solution. 

These coordinates have different definitions for t > r (we will refer to this inside 
the forward lightcone") and for t < r (that is, outside the lightcone). We have 

i: = rcosh?7, r = rsinh?7 (15) 

inside the lightcone and 

i: = rsinh?7, r = rcosh?7 (16) 

outside the lightcone. 

In the next section we use these coordinates in order to find new solutions of the hydro- 
dynamical equations. We investigate an ansatz for the flow and pressure field, which leads 
to a new and interesting class of exact, explicit and simple solutions. Among these solutions 
we find ones that overcome the shortcoming of almost all presently known exact solutions: 
the lack of acceleration. 

DERIVING NEW SOLUTIONS 

We parametrize the velocity with the f2(r, 77) rapidity as f = tanhf2(r, rj). In this section 
we take t > r, that is we are inside the forward lightcone. The Rindler variables r and 77 are 
defined by eq. (15). The Bjorken-solution is easily written down in terms of these variables 
saying Q = rj, and ^ = ^. Our new ansatz is the following generalization of the Bjorken 
flow: 

n = X7], A G M. (17) 

Here A is an (up to now) arbitrary constant, a parameter of the solution. Substituting (17) 
into the hydrodynamical equations (the Euler and energy conservation equations with the 
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e = up EoS) after some calculations detailed in Appendix A we obtain that p must have the 
following form: 

" = eo**-^'-)((A-l),)- 
with K being a constant. The value of A and K are constrained by the equation 

KK+{d- 1) ^^^^ cosh ((A - 1)7]) + 

sinh?7 

+ A (cosh^ ((A - l)ri) - k sinh^ ((A - = 0. (19) 

The detailed calculations which lead to these equations are found in Appendix A. 

Let us emphasize, that eq. (19) is derived from the hydrodynamical equations including 
the equations of state within lightcone, |t| > |r|, using Rindler coordinates. Outside the 
lightcone, for |t| < |r|, this equation is modified to eq. (69), that governs the "external" 
hydrodynamical solutions. 

In the next subsection we summarize and investigate the possible solutions. 



NEW EXACT SOLUTIONS 



There are 4 different sets of the parameters A, d, n and K which solve eq. (19), as discussed 
below. A fifth solution in this class is the A = 1 particular case, which is known as the Hwa- 
Bjorken solution in 1+1 dimensions, and whose multi-dimensional generalizations, often 
called as Buda-Lund type of relativistic solutions, were discovered in the past years [9, 10]. 

In Appendix B, we present the proof that these five solutions are the only possible non- 
trivial solutions in this class (where VL = Xt].) In Appendix B, we also prove the uniqueness 
of these solutions in a broader sense: If one allows A to be not only a constant but a proper- 
time dependent function (i.e. if f2 = A(r)?7), no other exact solutions exist, except the five 
solutions a.) - e.) listed in Table I. The properties of these solutions are detailed below. 



Flow profile and pressure 

In all cases, the velocity field and the pressure is expressed as 



V = tanhAr/, (20) 
P = P„(^)"'*'(cosh|)-"-"*\ (21) 
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Case 



a. j 

b. ) 

c. ) 

d. ) 

e. ) 



G 
G 
G 



d 
1 

4d-l 





ft 
K 






TABLE I: The new family of solutions. 

Table I shows the possible cases: every row represents one solution. In what follows, we 
discuss the various cases separately. 

• Case a.): In this case, A = 2, and the solution is obtained in arbitrary number of 
spatial dimensions, {d is arbitrary), but then k, should be equal to d. For instance, 
in three dimensions k should be 3, which is the EoS of an ultrarelativistic ideal gas 
or a photon gas. The corresponding flow profile has a non-vanishing relativistic ac- 
celeration: Uyd^u^ 7^ 0. We will see that this solution in the same form solves the 
hydrodynamical equations outside the lightcone, too. If we write down the velocity 
field in Minkowskian coordinates, it has the form both for t < r and t > r 

and we can solve the equation of motion for an individual particle, that is, we can 
determine the flow trajectories r{t). For a fluid element that is at position ro at to 
time-point, we have 

r(t) = -(v/l + (aot)2 + l) (23) 
ao 

These trajectories have constant ao acceleration in the local rest 



with ao 



2ro 



'o "-ol 



frame. A possible interesting application of this is mentioned briefly later. The tra- 
jectories of this solution are shown in Fig. 1. It has a realistic equation of state, 
corresponding to a massless ultra-relativistic gas in d dimensions. 

Note that we have found this solution flrst, using other methods. This solution can 
also be found using the criteria that the pressure depends only on proper-time and the 
solution should work in any d number of dimensions. These criteria select this case 
and the well known the A = 1 Hwa-Bjorken solution, detailed as case d.). 



• Case b.): In this case A = |, the number of spatial dimensions is arbitrary, but n = 1. 
The flow is accelerating. This case together with case c.) was found first in d = 3 
dimensions by T. S. Biro [12], and we generalized it for arbitrary d. An important 
point is that this case and case d.) are the only ones where the pressure field can have 
explicit ?7- dependence, and p can be finite in t] (actually, it is finite if and only if d 7^ 1). 

• Case c): This case is characterized with A = |, d is arbitrary, but k = ^^j^, that is, 
for example for (i = 1 we get k = 1, and for d = 3 it gives k = This case also was 
found in d = 3 by T. S. Biro first, and we have generalized it. The pressure field is 
also finite in 7] for d 1. 

• Case d.): In this case X = 1, d and k can be arbitrary. For d = 1 this case is known as 
the Hwa-Bjorken solution, and for d > 1 these solutions were found in refs. [9, 10]. The 
flow proflle in the d > 1 case is called Hubble flow. Thus these solutions were already 
known inside the light cone, but we quoted them for the sake of completeness. They 
are also the basis of new solutions, the the extensions of these solutions to the region 
outside the lightcone. These solutions are discussed in the next subsection. Inside the 
light-cone, this flow is not accelerating, u^d'^u'^ = 0. 

• Case e.): In this case, the value of the parameter A is arbitrary, but d and k must 
be equal to 1. The flow is accelerating, Uyd'^u^^ 7^ if A 7^ 1. The velocity fleld 
is pretty general, but the EoS is very particular. Nevertheless, this solution can be 
considered as a ,, smooth extrapolation" between the previous cases, although this kind 
of generalization — the arbitrary A case, which can be expected knowing the previous 
cases — works only with this special choice of d and k. In the forthcoming we shall 
use this solution to calculate the rapidity distribution. We show in Appendix C, that 
the width of the observable rapidity distribution is controlled by the parameter A. 
Hence in principle, the value of the parameter A of the hydrodynamical solution can 
be obtained from measurements. 

Outside the lightcone 

The hydrodynamical equations can also be rewritten in Rindler-coordinates outside the 
lightcone (that is, for \t\ < \r\). We investigate here how the solutions presented in the 
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FIG. 1: (Color online) Fluid trajectories of the A = 2 solution. 



previous subsection can be extended to this domain. The equations are found in Appendix 
A, here we only go through the possibilities. The statements below can be directly checked 
if one substitutes these expressions into eq. (69) in Appendix A, which is derived from the 
hydrodynamical equations in the \t\ < \r\ region. 

• Case a.): This solution can be extended outside the lightcone in a straightforward 
manner. The properties are the same: it works in arbitrary d with n = d, the velocity 
field has the same expression in Minkowskian coordinates as inside the lightcone (see 
eq. (22)), hence it is also uniformly accelerating. Fig. 1 shows the trajectories outside 
the lightcone, too. These trajectories are given by the same form, eq. (23), outside as 
well as inside the lightcone. 

• Case b.) and c): These solutions can be extended to t < r only in the d = 1 case, but 
these extensions are special cases of the extension of case e.) discussed below. These 
extensions do not depend explicitly on rj. (Such a dependence was only present for 
t > r only if c? 7^ 1.) 

• Case d.): An interesting new result is that there are solutions outside the lightcone 
which can be viewed as the extensions of the A = 1, or Hwa-Bjorken internal solutions. 
The velocity and the pressure have the form 



V = tanh?7 = -, 
r 



(24) 




(25) 
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Moreover, this extension is an accelerating flow. The equation of the trajectories is 



for ro > to- These trajectories are also of constant Oq acceleration, as the a.) (A = 2) 
case. However, the expression of Oq is different from that of case a.), we have now 



It is interesting to note at this point the following. As both the A = 2 and the A = 1 
external solutions are uniformly accelerating (in local rest frame), they both contain 
event horizons. This property of the uniform acceleration was utilized recently by 
Kharzeev and Tuchin [13] to describe thermalization in heavy ion reactions via the 
Unruh effect. 

• Case e.): The case when A is arbitrary and d = k = 1 can be extended to t < r 
without difficulties, the outside solution is also an accelerating one. The velocity and 
the pressure is given by 



but remember, that the definition of r and t] are given by eq. (16), which differs from 
the case of \t\ > \r\. 

Temperatures and densities 

The previously presented pressure and velocity fields solve eqs. (5) and (6), the Euler and 
the energy conservation equations. It is specific to the e = up EoS utilized here, that the 
continuity equation can be treated separately. But if we investigate the temperature and 
perhaps other charge densities, the relation between e and p does not fully determine the 
EoS. In this subsection we solve the continuity equations for some cases. 

Before doing this, we introduce the scaling variable S (r, f]) with the definition that its 
comoving derivative vanishes, that is (for one-dimensional or spherical flows) 




(26) 




V = tanhAr/, 



(27) 
(28) 



d^ dS dS 



(29) 



dt dt dr 
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S (r, 7]) 


A/1 


A = 1 


t> r 


(ro/r)^-isinh((A-l)77) 


7] 


t < r 


(ro/r)^~^cosh((A - l)r/) 


t/tq 



TABLE II: The S (r, rj) scaling functions inside and outside the hght-cone. Note, that the A = 1 
case has to be handled separately from the other cases. 



Every specific flow profile has a corresponding S function, which can be determined by 
solving eq. (29). Table II contains the expressions of S for the investigated class of solutions. 
(A remark: S is not uniquely determined by eq. (29), if 5* is a good scaling function, then 
any function of S is also a good scaling function. But we shall see, that this ambiguity does 
not really matter, as it cancels from the observables.) 

After this preparation, let us describe the solutions to the entropy (or in general, conti- 
nuity type of) equations. 

First, if there are no charges present at all, then from the (9) Gibbs-Duhem relation we 
have Ta = (k + 1) p, and it can be verified, that the solution of the entropy conservation , 
eq. (11) is 



o- = (To ( — ) J^aiS), 
\Po/ 



T = To 



P 



(30) 
(31) 



where p (r, rj) is given by eq. (21), and of course (k + 1) po = Tqcto. The function UaiS) > 
can be chosen arbitrarily. This approximation (that is, the baryonic charges are negligibly 
small) is often used for the description of high energy heavy ion collisions and of the quark- 
hadron transition of the early Universe. 

Second, if we have one conserved charge n, and as an ideal gas EoS we have p = nT, then 
the solution of the conservation equation is very similar: 



n = riQ 
T = To 



1 



(32) 
(33) 



The constants are chosen so that po = tiqTq. 
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The case of more than one charges can be treated similarly: we have {n + l)p = Ta + 
"^inifii, and the solution is 



K + l 



a = aol^] u^S), (34) 
\Po/ 

T = To (^) (35) 

\PoJ J^aiS) 

ni = nf u,{S), (36) 

\PoJ 

,W / P 



/i, = ^ M,{S), (37) 

\PoJ 

where the functions and the constants are constrained as Vi{S)M.i{S) = 1 and (k + l)po = 
^oC"o + X]i '^o Vo''- These forms are solutions to the hydrodynamical and continuity equations 
simultaneously. 

Let us emphasize, that the scaling functions of the scaling variable S are quite arbitrary, 
they are only constrained by using a physical normalization, which implies i^iiO) = M.i{0) = 
1. As we mentioned earlier, any function S' of the scaling variable 5* is again a good scaling 
variable, as long as 5"(0) = 0, because the same solutions are obtained when new scaling 
functions are chosen, that correspond to the new scaling variable as uKS'^S)) = Vi{S). Thus 
it is possible to use the most natural form of the scaling variable S, when describing these 
solutions. 



THE RAPIDITY DISTRIBUTION 

If we want to apply our solutions to the description of some physical phenomena, we have 
to calculate the final state observables from these solutions. We have chosen the rapidity 
distribution, and in order to be able to fit measured data, we apply here our new solution 
described as case e.) in the previous section. Here the parameter A (which we can call 
,, acceleration parameter", because it somehow influences the acceleration of the flow) can 
be arbitrary, but the dimension is restricted to = 1, and what is more inconvenient, k is also 
equal to 1. But the A parameter can be fitted to measured datasets, and this compensates 
the drawback of the specific EoS. 

The detailed calculation of the rapidity distribution is presented in Appendix C, here we 
briefly summarize the main steps. We use Boltzmann approximation, and neglect baryonic 
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charges and chemical potentials, so T and a are given by eqs. (30) and (31). Although 
the solution is valid only for the EoS of massless particles, we assume that at the freeze- 
out particles of rest mass m appear (e.g. pions with m = 140MeV). We have chosen the 
freeze-out condition as follows: the freeze-out hypersurface is pseudo-orthogonal to the four- 
velocity field u'^, and the temperature at 77 = reaches a given Tf value. The equation of 
this hypersurface is 

[^y~^ cosh {{X- 1)7]) = 1, (38) 

where Tf is the proper-time coordinate when the temperature reaches Tj at 77 = (see 
Appendix C for details). 

With a saddle-point integration in rj, for A > 1/2, m/Tf ^ 1 and V(j{s) = 1 we got 



cosh^2~^ - e ^/^ \ (39) 

y=0 \aJ 



dn dra 
d|/ d|/ 

with a = ^zY"- In an actual fit to data, the parameter a can be determined more directly, 
than A. From the value of the shape parameter a, the parameter of the acceleration A can 
thus be determined from the inverse formula as 

A = (40) 

This equation also indicates, that the fiat rapidity distribution corresponds to A — > 1 and 
a 00. 

In the formula eq. (39) and in the subsequent discussion of this distribution the upper 
sign always means the true 1 + 1 dimensional case (that is, not only the flow, but the 
produced particles have only one degree of freedom), while the lower sign is for the case 
when the 1 + 1 dimensional solution is embedded in the 1+3 dimensional space-time, and 
the produced particles have transverse degrees of freedom as well. The embedding means 
that the fluid occupies infinite volume, and it is homogeneous (it has no velocity, no pressure 
gradient components) in transverse directions. 

The Gaussian ,, width" of the eq. eq. (39) distribution (obtained from the second deriva- 
tive at y = 0) is 

Ay2 = , ui) 

^ m/Tj + l/2 + l/a ^ ^ 

We see from the expression of Ay"^ that it changes sign if A (and thus a) varies. It has two 
roots in A: Ay^ = if A = 1, or A = | ^1 + 2rn^fpfJ^ ■ If = 1 (this is the Hwa-Bjorken 
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FIG. 2: (Color online) Normalized rapidity distributions from the new solutions in 1+1 dimensions 
for various A, Tj and m values. Thick lines show the result of numerical integration, thin lines the 
analytic approximation from eq. (39). For A > 1 and not too big Tf it can be used within about 
10 % error. 

case), the distribution is entirely flat, and if A equals the other root, it is not flat but has only 
higher order non-vanishing derivatives at y = 0. In other cases, the sign of Ay"^ determines 
the behavior of the rapidity distribution (around y = 0): it has a minimum at y = 0, if 
Ay^ < 0, and it has a maximum (that is, it is approximately Gaussian) if Ay^ > 0. The 
typical cases are plotted in Fig. 2. This figure is just an illustration, the parameter values 
on this figure are not realistic ones, they were chosen in order to make the behavior of the 
distribution transparent in the different cases. 



APPLICATIONS 



In the previous section we have calculated the rapidity distribution from our new solution. 
We are now able to fit real rapidity measurements with our calculation, and we can extract 
the A parameter from the data. (Fig. 2 shows that the analytic approximation for not too 



16 



high temperatures has less than 10% error, so it can be used rehably.) We can also see on 
Fig. 2 that in the A — > 1 limiting case the distribution approximates the Bjorken-like flat 
one, and the more A differs from 1, the more the distribution deviates from the flat shape. 
That is, the acceleration effects make the distribution a finite one, and since our calculation 
takes acceleration into account, it provides a realistic description of the real distributions, 
which have finite widths. We present two applications: the improvement of the Bjorken 
energy density estimation, and an advanced reaction life-time estimation. 



Energy density estimation 

We follow Bjorken's method [5] and deviate from it at the point when acceleration effects 
(which are neglected by the Bjorken estimate) become important. Let us focus on a thin 
transverse piece of the produced matter at mid-rapidity, illustrated by Fig. 2 of ref. [5]. 
The radius R of this slab is estimated by the radius of the colliding hadrons or nuclei: 
R = l.lSA^/^fm. Its volume is dV = (-R^7r)rd?7, where r is the proper time of observation 
and drj is the space-time rapidity element corresponding to the slab. The energy content 
in this slab is dE = {■mt)dn, where (mj) is the average transverse mass at mid-rapidity, so 
similarly to Bjorken, the initial energy density is 

_ {rrit) dn 

Here tq is the proper-time of thermalization. (Bjorken's estimate was tq ^ Ifm.) For accel- 
erationless, boost-invariant Hwa-Bjorken flows Vo = Vf = Ui however, for our accelerating 

dy drif 
drjf drjo 

— (O 

dm 

estimation: the flrst one is the effect of work done by the pressure, and the second one 
characterizes how acceleration influences the expansion of the initial volume element. For 
an accelerating flow both factors should be greater than 1, which is indeed the case if A > 1. 
Thus the initial energy density can be accessed by an advanced estimation e^. as 

i^ = (2A-l)(^)"'. .« = J^^. (43) 

EBj \ToJ [R^TTjTody 

Here Sbj is the Bjorken estimation, which is recovered if ^ is flat (i.e. A = 1), but if A > 1, 
Eq is under- estimated by the Bjorken formula. Fig. 4 shows our flts to BRAHMS dn/dy 
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solution we have to apply a correction factor of -p^jr^. From our A G M solutions the shift 
of the point of maximum emittiviy is ^ = (2A — 1), while the volume element change is 
— — (rj/ro)^~^. These two factors contain the acceleration effects on the energy density 



FIG. 3: (Color online) Fluid trajectories of the A = 1.2 solution, which fits to BRAHMS ^ 
data. We illustrated the possible hypersurfaces (the hypersurfaces pseudo-orthogonal to the flow 
trajectories), among these ones there is the proper freeze-out hypersurface, selected by the criterion 

r(r? = o) = r^. 

data [14]. From these fits we have found A = 1.18 ± 0.01. On Fig. 3 we plotted the fluid 
trajectories and illustrated the ensemble of possible freeze-out hypersurfaces, as given by 
eq. (38) inside the lightcone for this specific A value. For illustration we also calculated the 
rapidity distribution for A = 1.18: the bottom panel on Fig. 2 shows it with this A value for 
pions, kaons and protons. 

Using the Bjorken estimate of Ebj = 5 GeV/fm^ as given in ref. [15], and Tf/ro = 8 ± 2 
fm/c, we find an initial energy density of Ec = (2.0 ± O.l)^^^ = 10.0 ± 0.5 GeV/fm'^. If the 
evolution deviates from a 1 + 1 dimensional perfect flow, then our estimation is only a lower 
limit for the initial energy density. 

Life-time estimation 

For a Hwa-Bjorken type of accelerationless, coasting longitudinal flow, Sinyukov and 
Makhlin [16] determined the longitudinal length of homogeneity as 



Here rrit is the transverse mass and tbj is the (Bjorken) freeze-out time. This result provides 
a means to determine the life-time of the reaction: it can be simply identified with tbj- 
However, if the flow is accelerating, the estimated origin of the trajectories is shifted back in 




(44) 
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a 


7.4 ± 0.13 


An 1 

Hi, V=0 
ay y 


294 ± 1 


x7ndf 


30.6/14 


CL 


0.6 % 


T/ (MeV) 


200 (fixed) 


m (MeV) 


140 (fixed) 


A 


1.18 ± 0.01 (derived) 



TABLE III: The parameters of the best fit with eq. (39) (1+3 dimensional case, lower signs) to 
dn/dy data of negative pions, as measured by the BRAHMS collaboration [14] in central (0-5%) 
Au+Au collisions at -^/s^vtv = 200 GeV. The free fit parameters were a and ^|y=o • The fit 
range was — 3 < y < 3 . The confidence level of the fit, CL was calculated from /NDF . The 
temperature parameter Tf was fixed to 200 MeV, corresponding to the slope parameter of the single 
particle spectra, as in the fitted 1+3 dimensional case a 1+1 dimensional exact hydro solution was 
embedded to a 1+3 dimensional space. For the pion mass, a fixed 140 MeV value was used. The 
parameter of the acceleration, A was calculated from the fitted parameter a using eq. (40). 




FIG. 4: (Color online) Panel (a): dn/dy data of negative pions, as measured by the BRAHMS 
collaboration [14] in central (0-5%) Au+Au collisions at y/s^N = 200 GeV, fitted with eq. (39) 
(1+3 dimensional case, lower signs). The fit range was —3 < y < 3, to exclude target and projectile 
rapidity region, CL = 0.6 %. The parameters of the fit are summarized in Table III. Panel (b): 
Ec/eBj ratio as a function of tj/tq. 
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proper-time, so the life-time of the reaction is under-estimated by tbj- (This was pointed out 
also in refs. [17-20].) From our solutions we have (for a broad but finite rapidity distribution, 
where the saddle-point approximation is valid): 



So the new estimation of the life-time, Tc contains a A multiplication factor. As mentioned 
in the previous subsection, BRAHMS rapidity distributions in Fig. 4 yield A = 1.18 ± 0.01, 
so they imply a 18 ±1 % increase in the estimated life-time of the reaction. 

SUMMARY 

In this paper, we presented new simple, explicit analytic solutions to the relativistic hydro- 
dynamical equations. These new solutions are characterized by a non-vanishing relativistic 
acceleration, which is a significant development compared to the already known solutions. 
Many of these new solutions work not only inside the lightcone, but outside of it as well. 
As a by-product, we have demonstrated that the well-known Hwa-Bjorken solution also has 
such an extension, and this extension turns out to have uniformly accelerating trajectories. 

We have calculated the rapidity distribution of the particles from our new solutions 
numerically, and presented a simple and reliable analytic approximation as well. These 
distributions have finite widths, hence they can be used to fit real experimental rapidity 
distribution data. We have fitted BRAHMS data, and obtained a value of A = 1.18 ± 0.01 
for the acceleration parameter of the solution. 

The famous Bjorken estimate of the initial energy density could be refined, taking ac- 
celeration effects into account. Using the Bjorken estimate given in the BRAHMS White 
Paper [15] {ebj = 5 GeV/fm^) and Tf/ro = 8 ± 2 fm/c, we find an initial energy density 
of Ec = (2.0 ± 0.1)eBj = 10.0 ± 0.5 GeV/fm^, which is a significant increase. We have also 
given an improved estimation for the life-time of the reaction: a 18± 1 % increase is due to 
the presence of acceleration, based on our model of the initial accelerating period. 
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Note that during the process of finalizing this manuscript, we have encountered two inter- 
esting, new family of relativistic hydrodynamical solutions, in ref. [23] and also in ref . [24] , 
both of which look to be very interesting and deserve more detailed, further investigations. 

APPENDIX A: HYDRODYNAMICAL EQUATIONS IN RINDLER COORDI- 



In this Appendix we present the derivation of eqs. (18) and (19) in more detail. First 
we rewrite the hydrodynamical equations in terms of the Rindler variables, then investigate 
the assumption for the velocity field made in eq. (17). In this Appendix, as in the body of 
the paper, Q stands for the rapidity of the flow, v = tanhfi, and here we use the notation 
^ = ^ + w^, this appears on the r.h.s. of the Euler equation. We treat the two different 
domains of the variables — inside and outside the lightcone — carefully, because the Rindler 
coordinates have different definitions. We shall see that the Q = t] case is a particular one. 
(We refer to it as Hwa-Bjorken case, but emphasize that for d > 1 such solutions were found 
only recently.) So we investigate it in a separate subsection, in the following two subsections 
— where we derive the general equations inside and outside the lightcone — we assume 



Inside the forward lightcone 

First, we deal with the t > r, t > case. The definition of the Rindler coordinates, the 
inverse formulas and the derivation rules are as follows: 



NATES 



A^l. 



r = — , 



t = T cosh r] , 



r = T sinh 77, 




(46) 
(47) 
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d d sinh?7 d , . 

d d cosh?7 d , , 

— = -sinh?7— H — , (49) 

or OT T or] 

d cosh [VL — rj) d sinh [VL — rj) d . . 

dt coshf2 dr rcoshfi 9//' 

d sinh {Q — rj) d ^ cosh (f2 — 77) d 
dr coshfi dr rcoshfi drf 

The domain of the variables is —00 < rj < +00, and < r < 00. We introduce the notation 

Q = -j^lnp. (The convenience of this notation is due to the special EoS we utilize here, 

the e = up one.) Rewriting and rearranging the Euler and energy conservation equations 

we obtain the following: 

on IdQ , , fdQ ldn\ ^ . . 

_ + __^ + tanh(n-„)(^ + --)=0, (52) 

dQ Idn d-1 sinh il 
dr T drj cosh (f2 — 77) r sinh rj 

+ tanh(n-,)(^ + ^^)=0. (53) 

Here we used the expression for the divergence of the velocity field: (Vv) = + ^^v. 

Our special assumption was that Q = Xr], so substituting ^ from eq. (52) into (53) we 
have 

tanh((A-l),)(r^ + A)+^ = 0, (54) 

.r^ + (rf-l)^cosh((A-l),) + 
OT smh rj 

+ A (cosh^ ((A - l)ri) - k sinh^ ((A - l)r])) = 0. (55) 

From eq. (55) we see that must be r-independent, that is Q must have the form 
Q = H (rj) + K (rj) In r. Then from eq. (54) we have that K (77) is also constant. Again with 
eq. (55) we have then 

Q = K\nT + H{r]), (56) 

where the H (rj) function is determined by eq. (54). The solution is easily obtained as 

K -\- X 

H iv) = — T 1^ cosh ((A — 1)77) + const. (57) 

A — 1 

for A 7^ 1 , so p has the form 

P = Po[-) — . K+x,^,, — • (58) 



r / 



cosh^('^+^) ((A - 1)7/) 
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The integration constant was absorbed into the definitions of the po and tq initial values. 
For the K constant we have from eq. (55) 

KK+{d- 1) ^^^^ cosh ((A - 1)7]) + 

sinh?7 

+ A (cosh^ ((A - 1)7]) - K sinh^ ((A - 1)7])) = 0. (59) 

Essentially this is eq. (19), so what we need is to find such A and K constants, which solve 
this equation, then we have the pressure and the velocity field. The solutions were found, 
they are presented in the body of the paper. 



Outside the forward lightcone 

In this case the calculation goes almost in the same way as inside the lightcone. The 
definitions and the derivation rules are: 

t = rsinh?7 , r = rcosh?7, (60) 

r = ^7-2 _ ^2 ^ ^ — arctanh-, (61) 

r 

d d cosh?7 d , . 

- = -smh,- + ^-, (62) 

d d sinh?7 d . 

— = cosh?7- — , (63) 

or OT T 07] 



d sinh {VL — 7]) d cosh (f2 — 77) d 

dt coshf2 dr rcoshfi d7]^ 

d cosh [VL — 7]) d sinh (f2 — 77) d 

dr coshf2 dr rcoshfi dr] 

We have in this case for the hydrodynamical equations 



(64) 
(65) 



dQ Idn , , fdn ldQ\ 

dVt ^ K dQ d — 1 sinh 



dr T dr] cosh {Vt — r]) t cosh r] 

+ tanh {n-r])(K^ + -^\=Q. (67) 



dr T dr] 

If we now assume ^2 = Ar/, then going through the same steps as inside the lightcone we 
have (for A 7^ 1) 



P=Po 



(^) ^^" sinh-^^'^^^) ((A - l)r^) , (68) 
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and the equation which constrains K and A is 



k{K + X) , , , sinh Xri 
(^-1) 



sinh ((A — 1)77) cosh?7 
+ (/€- 1) Asinh((A- = 0. (69) 

This equation is analogous to eq. (59), but it is shghtly different. We used this equation to 
investigate the solutions outside the lightcone. 

The Hwa-Bjorken case 

The A = 1 case is the Hwa-Bjorken solution in 1+1 dimensions, and the Buda-Lund case 
(refs. [9, 10]) in 1+d dimensions, which is simply obtained from eqs. (52) and (53). We have 

dQ ^ dQ d 



^ = V, P = Po(-) . (71) 



and the solution is 

f2 = 77. V = Vn ( ^ 

What is perhaps more interesting, is that this solution can be extended to outside the 
lightcone, as we see it from eqs. (66) and (67), if we substitute A = 1: 

|3,i.O, f.y-^u.,,, (72) 

or r orj K 

so the ,, outside Hwa-Bjorken-Buda-Lund solution" is 

n = r/, p = po [-^) cosh- — (''+^) rj. (73) 
This is an accelerating solution, its properties are discussed in the body of the paper. 

APPENDIX B: TAYLOR EXPANSION AND UNIQUENESS 

In this Appendix we sketch the proof of uniqueness of the presented solutions. We deal 
only with the t > r case. We investigate a less restrictive assumption for the velocity field 
than in the body of the paper, namely allow the A factor to be proper-time dependent: 

n = X{T) 7]. (74) 
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We denote the derivation operator by a comma, and if possible, we suppress the r 
argument of the functions in the notation. The A(r) = 1 case was investigated in Appendix 
A as well as in the body of the paper, so we assume A(r) 7^ 1 throughout this Appendix. 
We will use a Taylor-expansion around rj = for Q and the pressure (that is, for Q = 

Q = Air)+Bir)v + ^V^ + ^v' + ^V^ + ..., (75) 

We have nothing more to do than to put this into eqs. (52) and (53), and make the l.h.s. 
vanish order by order in powers of 77. It can be done without any essential difficulty, but it 
is a very long calculation. Here in this Appendix we only consider the case when one makes 
an additional assumption for the pressure field, namely, that it factorizes in terms of r and 
r], as in the solutions presented in the body of the paper: 

P = F{t)H [t]) (76) 

That is, for the Taylor expansion we have 

Q = AiT) + Br^ + ^v' + jv' + jv' + ---, (77) 

where A (r) = -j;^ In F (r) and the other functions are constants (denoted by the same 
letter), the expansion coefficients of \nH (77). We will need only the second order terms 
for the present calculation: we expand eq. (52) to 0{r]^), and eq. (53) to (9 (77^). The 
following auxiliary formulas are needed: 

tanh ((A - 1)7]) = {X-l)r]-^{X-l)'r]^ + ..., 

sinh Xr] 1 ^ 1 ^ ^ ^ t \ r\ ON 2 



A - -A (A -1) (A -2) 7/' 



cosh ((A — 1)77) sinh 7/ 3 
+^A (A - 1) (A - 2) (6A' - 12A + 7) r/' + . . . 



B + {X' + C + {X - 1) {A' + X)}r] + Dr]^ + 



Eq. (52) thus gives 



+ \^E--{X-lYiA' + X)jr]^ = 0, (78) 

and eq. (53) is (to second order in rj) 

{kA' + dX} + (A - 1) kBt] + 
+ (A-1){^^A(2-A) + A' + kC}7/2 = 0. (79) 
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We immediately see that B = D = 0, and we obtain the following equations: 

KA' + d\ = 0, (80) 

A' + C + (A-1)(A' + A) = 0, (81) 

(2 - A) + A' + = 0, (82) 

3 

^ (A - 1)' iA' + X)-E = 0. (83) 

For brevity we did not write the O (77^) term in eq. (79), but we will use the equation, which 
follows from the condition that the O (77^) term vanish: 

^A (A - 1) (A - 2) (6A2 - 12A + 7) + 

+ {X-1)kE- {kC + A') = 0. (84) 

Now substituting A' from eq. (80) into eq. (81), and then A' from eq. (81) in eq. (82), this 
last equation will contain only A (not its derivative), and the C constant. So eq. (82) will 
yield an algebraic equation for A, that is, A must be a constant, thus A' vanishes. It follows 
that A' is also constant. So the constants are 

A' = --A, (85) 
C=-A(A-l)(l-^), (86) 

E^ix(^-lf{l-^}), (87) 

and from eqs. (82) and (84) we have two equations containing only A: 

2 (A - 1)' {K-d) + (A - 2) (6A2 - 12A + 7) = 0. (89) 

We excluded the trivial A = case, and as mentioned before, the A = 1 Hwa-Bjorken case 
was treated separately. Substituting eq. (88) into eq. (89) we obtain 

{d - 1) (A - 2) (4A2 - 8A + 3) = 0. (90) 

It is now easy to see that there are only those solutions which we presented in the body 
of the paper: if ci = 1, then eq. (88) implies that k = 1, but then eq. (90) is satisfied by 
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any value of the parameter A. Second, if A = 2, then eq. (90) is satisfied in any number 
of dimensions but eq. (88) imphes a condition for the equation of state, n = d. Third, 



other solutions of eq. (88). So the only possibilities are those listed in Table I. They are 
indeed solutions, as can be verified straigthforwardly by substituting these solutions to the 
hydrodynamical equations in Rindler coordinates. The conclusion of this Appendix is that 
we cannot extend the investigated class of solutions in an easy way even if we allow A to be 
r-dependent. 

APPENDIX C: RAPIDITY DISTRIBUTIONS 

In this Appendix we derive eq. (39), the formula which gives the rapidity distribution of 
the produced particles for our new 1+1 dimensional solutions, which are given by eqs. (20) 
and (21), and the constants are those in case e.) in Table I. We assumed = 0, that is 
e + p = Ta. We utilized the following expressions for cr and T, similarly to eqs. (30) and 



so we assumed i^aiS) = 1 for the scaling function, and p(r, rj) is given by eq. (21). Of course, 
Pf = cTfTf. The subscript / means freeze-out. (The constants in eqs. (30) and (31) can 
always be re-scaled, so we denoted the proper-time coordinate of the freeze-out hypersurface 
at ?7 = by Tj, and Pf,Tf and aj are the values of the thermodynamical quantities in this 
space-time point.) Here, as in the body of the paper, m is the mass of the produced particle, 
and we parametrize the momentum of it as k'^ = {E,p) = m (coshy, sinhy), with y being 
the rapidity. 

We calculate the rapidity distribution in two cases. First, in the fully 1 + 1 dimensional 
case, when not only the flow but also the produced particles have only 1 degrees of freedom, 
as they propagate along the z axis. Second, in the case, when this 1 + 1 dimensional solution 
is embedded in a 1 + 3 dimensional space, and the produced particles can have transverse 
momentum as well. 



the other two roots of eq. (90) are A = |, and A = |, and it is clear, that there are no 



(31): 




(91) 
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The 1 + 1 dimensional case 



The rapidity distribution is given by 

^ = E— = / d2x5(a;^F), (92) 
dy dp J 

where S{x^, k'^) is the source function, and the notation d^x refers to the fact that we are 
in an 1 + 1 dimensional space-time. The exphcit form of S'(x^, k'^) is 

^(^- - in iw) - ^) - 

where the notations are: u'^ is the four-velocity of the flow, fi is the chemical potential, g 
is the spin-degeneracy factor, and dS^(x) is the Cooper- Frye flux term. (The Cooper- Frye 
flux term, together with the delta function, which selects the freeze-out hypersurface, is 
essentially the Lebesgue vector-measure of the freeze-out hypersurface). 

In our particular solution, /i(x) = 0, and k^u^{x) = m cosh {Xt] — y) . We selected the 
freeze-out condition in a way to take advantage of the simplicity of this relation: our con- 
dition is that the temperature at ?7 = reaches a given Ty value, and the hyper-surface is 
pseudo-orthogonal to the four-velocity field, or in other words, u^{x) is parallel to dS'^(x). 
(This hypersurface is not of constant temperature nor of constant r. This choice is also moti- 
vated by the success of the Buda-Lund model in describing particle spectra, see e.g. [21, 22].) 
Another formulation of this condition is that the freeze-out hypersurface is ,, locally synchro- 
nized" : an observer which moves together with the flow sees this hypersurface as a constant 
time surface in its neighborhood. The equation of this hypersurface is 

(^)^~^cosh((A-l)r/) = 1. (94) 

We can now calculate the Cooper-Frye term. We have for the integration 

d^x5 {x - Xf) k^d'L^'ix) 

/oo ^ 
dr] rriTf cosh (A?7 — y) cosh>^~^ ((A — 1)?]) . (95) 
-oo 

So the rapidity distribution is 

dfl /"OO jjiTf 1 

— = / d?77-— 7 cosh (A?7 — ?/) cosh^^~ ((A — 1)?]) x 
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X exp I — — cosh (A?7 — y) cosh>^ ((A — l)ri) I . (96) 
V / 

This integration can be calculated numerically, but an analytic approximation is also avail- 
able. The formula of the saddle-point integration is 



dxf{x)g{x)^g{xo)^^^j^, (97) 

for an f{x) function, that has a sharp maximum peak at xq, and a g{x), which changes 
smoothly around xq- Essentially this is an approximation where one treats the smooth func- 
tion as a constant g{xQ), and the sharp function as a Gaussian peak with parameters obtained 
from /(xq) and /"(xq). This method can be applied very effectively when the peaked func- 
tion has a parameter, which rules its half-width, so this parameter can be changed in such 
a way that f{x) approximates a Dirac-delta. The more f{x) is like a Dirac-delta, the better 
the saddle-point approximation is (and in the limiting case it becomes exact). 

In our case, /(?]) is the exponential, g{ri) is the remaining, and ^ is the parameter, which 
governs the half-width of /(?]). For convenience, we introduce the a = ^3^" notation. We 
have 



a 



f'iVo) = -fiVo)^X (2A - 1) cosh^-i (I) , (99) 

g{Vo) = '^cosh^^(y). (100) 
Putting all this together, we obtain (approximately) 

= A-co.ht-. (i) exp cosh" (f ) + ^} , (101) 

This is the expression of 4^ in the true 1 + 1 dimensional case. The normalization constant 



dn 



dy 

, the value of ^ at mid-rapidity. 



The 1 + 3 dimensional, embedded case 

If we take the transverse dynamics into account, that is, we embed the 1 + 1 di- 
mensional solution in a 3 dimensional space (with Vx = Vy = and no pressure gra- 
dients in transverse directions), the previous (101)-(102) formulas need to be modified. 

29 



We have (in 3 spatial dimensions) ^2nhf instead of and instead of m, where 
m-r = E'^ — -p^ = ^Jm^ + + transverse mass. With these modifications, the 

above (101)-(102) expressions have the following meaning: 

„ dn dn . ^ , « i / V \ 

E—— = = ATcosha - X 



d^k dy2'irmTdmT V« 

xexp|-!^co=h-m+!^l. (103) 



Tf \a/ Tf 



^'mTA(2A-l)(27r;i)' \ Tj ' ^ ' 

In order to obtain the rapidity distribution, we have to integrate over the transverse mass: 



dn ^ / 27rTf niTTf 

—— = 271 I drrLTrriTx , , / 

Jm mT\{2\-l){2Txhf 



xexpj-^cosh" ©• 



(105) 

Here m is the real mass of the produced particles, e.g. m = 140 MeV for pions. The integral 
over mx can be also calculated in the limiting case of ^ oo. (This is the same condition 
as that for the saddle-point integration to be exact.) In this case, we can apply the simple 
approximation 

h{x) exp {-(3x) dx ^ ^ (106) 

for a -I, and for any /i(x), that varies smoothly around a. Our case is h{mT) = m'^'^ . So 
we have 

^ = AT'cosh-t-i (I) exp 1-^ cosh" a) ] , (107) 

U' = X -— — 108 

V A(2A - 1) 2Txh' ^ ' 

with A being the transverse cross section of the fluid. These are the results quoted in 

eq. (39). 



APPENDIX D: GENERAL SOLUTION FOR THE d = 1, k = I CASE, INVESTI- 
GATION OF STABILITY 

If e = Kp, and we take the choice of k = 1, then the hydrodynamical equations (that is, 
eqs. (5) and (6)) have an intrinsic simplicity. Let us write them down in this case! In this 



30 



Appendix we denote the derivation with respect to t by dotting, and the derivation with 

respect to r by comma. So the Euler equation is 

1 1 
-—^ {v + vv') = -—{p' + vp) . (109) 

We have e + p = 2p for the enthalpy density. The energy conservation is 

^ (p + vp') = -^(^^ {{vy+v {v') ') - v', (110) 

or in another form 

^—^y' + vi,) = -}-^p + Vp'). (Ill) 

1 — f ^ 2p 

We used the fact that in d = 1 dimensions (Vv) = v' simply. Taking a proper linear 
combination of eqs. (109) and (111) we have the following: 

n' = -^p, (112) 



1-^2 2p 



\ — v'^ 2p' 

Here we again used the Q = arctanht; notation. These equations are analogous to the 
Cauchy-Riemann equations known in complex analysis, and the general solution of them 
is easily written down: it contains two arbitrary functions, F {t + r) and G {t — r), one for 
each direction of wave-propagation: 

n{t,r) = i{G(t-r) + F(t + r)}, (114) 
\n{p (t, r) /po) = {G{t~r)-F{t + r)} . (115) 

Thus this case is a particular one: the general solution of the relativistic hydrodynamical 
equations is found, and then it can be easily fitted to a given initial condition. If we prescribe 
p{t{r),r) and Q{t{r),r) on some t = t{r) Cauchy hypersurface, then the proper F and G 
functions can be found as detailed below, and then the solution of the initial value problem 
is that the two wave-shapes propagate in the opposite direction. As a first example, let us 
consider the Q = Xt] solution presented in the body of the paper. It can be easily verified, 
that the following F and G functions indeed yield it: 

(116) 
(117) 
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To simplify the presentation of the general case, let us introduce 

n(t,r) = ^ln(p(t,r)/po), (118) 

and give the initial boundary conditions in the form of three functions, the Cauchy hyper- 
surface to{r), the initial velocity field v = tanh(f2) given by the form of Q{to{r),r) and the 
initial pressure distribution is given hj p = poexp(2Il{to{r),r). 

Let us introduce the positive and the negative light-cone coordinates on the initial Cauchy 
hypersurface as 

x^{r) =to{r) ±r. (119) 

It is important to assume, that on the Cauchy hypersurface, these relations can be inverted, 
and r can be expressed as a function of both the positive and the negative light-cone coor- 
dinates: 

r = /(x+), (120) 
r = g{x~). (121) 

For example, if the initial Cauchy hypersurface is given by r = tq = const, then the equation 
of this hypersurface in time and coordinate is to{r) = ^/t^ + r'^ and the relation between 
the coordinate r and the light-cone coordinates, 



x^{r) =to{r)±r = + r^±r, (122) 
which relations can be inverted as 

+ 2 

r = = T - 2^' ^'2') 

r = g{x-) = + (124) 

Then the general solution of the hydrodynamical equations in this case simplifies in terms 
of Vt and 11 to the following relations: 

fi(t, r) = ^ [fi(to(/(t + r)), fit + r)) + fi(to(^?(t - r)),g{t - r))] + (125) 

+ ^ [n(to(/(t + r)), f{t + r)) - Ii{tMt - r)),g{t - r))] , (126) 

n(t, r) = ^ [-n{t,{f{t + r)), f{t + r)) + ^{tMt - r)),g{t - r))] - (127) 

-\ [n(to(/(t + r)), f{t + r)) + Ii{tMt - r)), g{t - r))] . (128) 
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A by-product of the treatment of this Appendix is, that any 1 + 1 dimensional relativistic 
hydrodynamical solution with e = p EoS (including our one presented in the body of the 
paper, with Q = Xr]) is always stable in terms of Q and 11: if there is a perturbation in the 
initial condition for IT and Q, then this perturbation propagates in both directions, but it 
does not grow exponentially. It is interesting to contrast this behaviour, the stability of our 
solutions in terms of rapidity and logarithmic pressure, to the numerically found instability 
of the famous Landau-Khalatnikov solution for the perturbation of their initial conditions. 
Such a stability of our new, general solutions makes them particularly attractive for testing 
various algorithms that integrate the solutions of relativistic hydrodynamics numerically. 
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